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^ Abstract: 

Usual analyses based on scans of the seesaw parameter-space can be biassed since they 
J> do not cover in a fair way the complete parameter-space. More precisely, we show that 

in the common "i?-parametrization" , many acceptable i?-matrices, compatible with the 
£ — . perturbativity of Yukawa couplings, are normally disregarded from the beginning, which 

produces biasses in the results. We give a straightforward procedure to scan the space 
of complex .R-matrices in a complete way, giving a very simple rule to incorporate the 
perturbativity requirement as a condition for the entries of the i?-matrix, something not 
considered before. As a relevant application of this, we show that the extended believe that 
BR(^ — > e, 7) in supersymmetric seesaw models depends strongly on the value of #13 is an 



"optical effect" produced by such biassed scans, and does not hold after a careful analytical 
and numerical study. When the complete scan is done, BR(^ — > e, 7) gets very insensitive 
to #13. Moreover, the values of the branching ratio are typically larger than those quoted in 
the literature, due to the large number of acceptable points in the parameter-space which 
were not considered before. Including (unflavoured) leptogenesis does not introduce any 
further dependence on #13, although decreases the typical value of BR(fi — > e, 7). 
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1 Introduction 

There is a generalized believe [1],[2] that, in the context of the supersymmetric (SUSY) 
seesaw scenario, the value of the neutrino mixing angle #13 has a strong impact on Lepton 
Flavour Violation (LFV) processes, in particular on the branching ratio BR(/u — > e, 7). 
The basic idea is the following. 

As is well known, starting with universal conditions for the soft mass matrices (and thus 
implementing minimal flavour violation), the renormalization group (RG) running induces 
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non- vanishing off-diagonal entries in the (left-handed) slepton mass matrix, hil 2 . Such 
entries are mainly generated by the matrix of neutrino Yukawa couplings, Y u ; more pre- 
cisely (rriL 2 )ij ~ (yJYj,)^-. Then these off-diagonal entries enable LFV processes through 
one-loop diagrams [3]. On the other hand, for given low-energy observables (i.e. neu- 
trino masses and neutrino mixing matrix), the matrix YiY u has plenty of freedom. This 
is because in the standard seesaw scenario there are more initial high-energy parameters 
(18) than low-energy observables (9). Consequently, for given low-energy observables the 
rate of LFV processes can vary within a certain range. The claim of ref. [2] is that this 
range typically shifts to larger values as the #13 mixing angle increases; furthermore, such 
behaviour is strengthened by leptogenesis constraints. The effect is apparently very im- 
portant, producing increases of several orders of magnitude in BR(^ — > e, 7) as #13 grows 
within its experimental window. This remarkable behaviour has been shown up mainly by 
performing scans upon the seesaw parameters (with and without leptogenesis constraints). 
An analytical explanation for it has been offered in [4]. 

In this paper we will show that these results are essentially an optical effect produced 
by scanning only a part of the whole parameter space of the seesaw. As we will see, when 
the whole parameter space is considered, the impact of #13 is very small, even negligible. We 
will also show that, typically, the scans of the seesaw parameters in the literature are both 
incomplete and biassed because they do not consider the whole freedom in the parameter 
space compatible with the requirement of perturbativity. Such lack of completeness is 
more important for some parameters than for others, thus the bias. We will present easy 
rules to perform complete explorations incorporating the perturbativity requirement. The 
mentioned common bias is the main reason for the observed conspicuous dependence on 



In section 2 we present the framework and set the notation. Section 3 is devoted to 
describe the two parameterizations of the seesaw we are dealing with in this work. In 
section 4 we show and explain the apparent contradiction between the two approaches, 
concerning the dependence of BR(/j — > e, 7) on the #13 angle. In section 5 we give a 
procedure to scan the seesaw parameter space in a fair way, incorporating the constraint 
of perturbativity (details are given in Appendix A). In the same section and section 6, we 
study both in an analytical and numerical way the dependence of BR(^, — > e,7) on 613, 
making use of the mentioned scan. Section 7 is devoted to the inclusion of leptogenesis 
constraints. Finally, in section 8 we present our main conclusions. 

2 Framework and notation 

From now on we use the conventions and notation of ref. [5]. 
In the standard SUSY seesaw the relevant superpotential is 



where L (e^) are the leptonic doublets (charged singlets), vr are the right-handed neutrinos 
and if 1,2 are the two supersymmetric Higgs doublets. Y e and Y u are Yukawa matrices 
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in flavor space (flavor indices are dropped) and M is the Majorana mass matrix of right- 
handed neutrinos. Once right-handed neutrinos are decoupled (at the seesaw scale ~ M) 
a mass operator is left in the effective superpotential, 

W eS D e% T Y e L-H 1 + ^(Y V LH 2 ) T M- 1 (Y V LH 2 ) . (2.2) 
Then, the effective mass matrix for the light (~ left-handed) neutrinos is 

M v = (H°) 2 k , (2.3) 

with 

k = Y V T M~ X Y V . (2.4) 

Note that the seesaw equations (2.2), (2.4) are valid at the seesaw scale. Besides, they are 
obtained using the (reasonable) approximation of decoupling at a unique threshold, instead 
of a (more accurate) decoupling in three steps (the three right-handed neutrino masses). 

The neutrino mass-eigenvalues, mi = {H 2 ) 2 Ki, and the mixing matrix, C/mnSj are given 

by 

At = t^MNS K UySNS, At = diag(«i,« 2 ,K3) ) (2.5) 
with K\ < K2 < K3. The standard parametrization of the MNS matrix is 
/ C13C12 C13S12 s 13 e~ l5 \ ( e' 1 ^ 2 



-C23S12 - S 2 3Si 3 Ci 2 e' lS C23C12 - S23Sl3Si2e l5 S23C13 
\ S23S12 - C23Sl3Cl2e lS -S23C12 - C 2 3Sl3Sl2e l<5 C23C13 



V 



e -i4>'/2 I (2.6) 



From now on, for the sake of notation clarity, we will drop labels in the neutrino mixing 
and Yukawa matrices: U = J7mns> Y = Y v . 



3 Parameterizations of the Seesaw 

The high-energy seesaw Lagrangian, given by the superpotential (2.1), is determined by 
the entries of the Y and M matrices, which contain 18 independent parameters. On 
the other hand, there are 9 low-energy neutrino observables: the three neutrino masses, 
oc K{, and the three mixing angles and the three phases contained in U. Hence, for given 
values of the low-energy observables, the freedom in the seesaw Lagrangian expands a 9- 
dimensional parameter space. There are two main ways of describing such space (or, in 
other words, of parametrizing our ignorance). We will call them the R— parametrization 
and the Vl— parametrization. 
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3.1 R— parametrization 



It was shown in [5] that, for a given set of low-energy observables, k% and U, the neutrino 
Yukawa matrix (at the seesaw scale) has the form: 



(3.1) 



where D ^jj = diagjv^MTj-, = diag{y / Ki}, and R is a complex orthogonal (3 x 3) 

matrix. So, the 9 see-saw parameters that parametrize our ignorance are the three right- 
handed masses and the 3 complex angles defining R. A usual parametrization of R is 



R 



( c 2 c 3 
C2S3 



-C1S3 - S1S2C3 
C1C3 - S1S2S3 
S1C2 



S1S3 - C1S2C3 
-S1C3 - C1S2S3 
C\C2 



(3.2) 



where (q) are the sine (cosine) of the three complex angles Eq.(3.2) is general up to 
reflections changing the sign of det R. 



3.2 Vl— parametrization 

An alternative to the R— parametrization is the so-called Vl— parametrization, see ref. [6]. 
For a given set of low-energy observables, Ki and U, the neutrino Yukawa matrix (at the 
seesaw scale) can be written as 

Y = VrDyVI . (3.3) 

Here Dy is the diagonal matrix containing the three (real and positive) neutrino Yukawa 
couplings, in; Vl is a unitary matrix with identical structure as the MNS matrix [eq.(2.6)], 
but, of course, with three different mixing angles and three different phases; and Vr has also 
identical structure but with the diagonal matrix of phases acting from the left (for more 
details see e.g. [7]). Here, the 9 independent parameters that parametrize our ignorance 
are the three Yukawa couplings, yi, and the three angles and three phases contained in Vl- 
The Vr— matrix and the three right-handed neutrino masses are obtained by substituting 
eq.(3.3) in the seesaw expression (2.4), namely 

V R D M V* R = DyVlUD~ l U T V£D Y , (3.4) 

i.e. Vr is the unitary matrix that diagonalizes the symmetric matrix in the right hand 
side of (3.4) and Dm is the corresponding diagonal matrix, which contains the three right- 
handed masses. 



4 The impact of #13 

As mentioned in the introduction, even starting with diagonal and universal soft masses, 
the RG running generates off-diagonal entries in them. In particular, the (left) slepton mass 
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matrix gets off-diagonal entries (mL 2 )ij, i 7^ j, proportional (in the leading-log approxima- 
tion) to the (YtY)jj matrix element. On the other hand, at first order in the mass-insertion 
expansion, the branching ratio of the LFV process Zj — > L-,7 (with Zjj charged leptons of 
the i,j families) is non- vanishing and proportional to the squared of the corresponding 
off-diagonal entry, 

Consequently, the dependence of, say BR(// — > e, 7), on the MNS matrix (and in 
particular on 9i 3 ) occurs mainly via the dependence of (Y^Y^i on it. This becomes more 
clear from the following approximate formula for the branching ratio, which will be useful 
later for qualitative discussions, 



(m L 2 )d 2 - 



BR(/i — > e, 7) 



Q 



G 2 F m% 



(Y t Y ; 



21 



tan 2 . 



(4.1) 



Here represents a typical supersymmetric leptonic mass, and mo, Aq are the universal 
scalar mass and the universal trilinear coupling at the unification scale Mx- 

Next, we analyze the dependence of (Y^Y^i on #13 using the two parametrizations of 
the seesaw discussed in sect. 3. We will find first a kind of "paradox", then we will discuss 
its explanation. 



4.1 A "paradox" 

In the R— parametrization the Y^Y matrix can be easily obtained from eq.(3.1): 

YtY = UD^D M RD^ . (4.2) 

This equation tells us that, for given right-handed masses {Dm) and a given R— matrix, 
Y^Y (and thus LFV processes) has a non-trivial dependence on U . Let us concentrate 
for the moment on the (Y^Y^i matrix element, which is the relevant one for \i — >■ e, 7. 
Assuming a hierarchical spectrum of neutrinos, k\ <^ K2 -C K3, we can expand (Y"l"Y)2i in 
powers of ^J~Ki. The first term of such expansion is 

(YtY) 21 = k 3 U 23 [r^D m r] 33 U1 3 + ••• (4.3) 



Since \U\ 3 \ = \si 3 \ and \U2 3 \ ~ l/v2, we see that, in this approximation, 

(YtY) 21 oc s 13 . (4.4) 

This dependence on si 3 is quite strong, and is really the source of the dependence of 
BK(fi — > e, 7) on ~ s 2 3 observed in the literature (this was also noticed in [4]). An 
analogous argument shows that BR(r — > e, 7) has a similar dependence on si 3 , while 
BR(r — > fi'y) is almost independent of #13. We will discuss soon the validity of the previous 
expansion, and thus of these results. 

Let us now turn to the Vl— parametrization. From (3.3) 

Y f Y = V L D\vl . (4.5) 
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Clearly, now Y'Y does not depend at all on U. For given Yukawa couplings (Dy) any 
choice of Vl is compatible with any choice of U and thus of #13. So, varying #13 does not 
affect YtY at all. This result seems to be in contradiction with the one obtained using the 
R— parametrization. 

One might argue that changing #13 in the Vl —parametrization means moving along 
a line of constant yi and Vl in the seesaw parameter space; while changing #13 in the 
R— parametrization means moving along a line of constant Mj and R. It may happen that 
(Y'Y)2i keeps constant along the first line but it increases along the second one. This is 
true, but even assuming this possibility there remains a conflict, as we are about to see. 

Let us work first in the Vl— parametrization. Assuming hierarchical neutrino Yukawa 
couplings, y\ <C y\ <C yf, it is obvious that the choice of Vl that maximizes (YTY)2i in 
eq.(4.5) is |(Vt)i 3 | = |04) 23 | = Vv^- Then 

/ \ max I 

(Y«Y) 2i = -„?. (4.6) 

As mentioned above, this value is available for any choice of U. Let us keep fixed all 
mixings and phases in U, except #13. For each value of #13, the Y'Y matrix remains 
the same, but Dm and Vr change according to eq.(3.4). The corresponding Y matrix is 
given by (3.3). Suppose we choose #13 = 3° and then calculate Dm and Vr, and write 
Y. This matrix can be easily expressed in the R— parametrization. Namely, in eq.(3.1) we 
can straightforwardly solve R in terms of Y, U, Dm and D K , which are known. We can 
wonder now what happens if we keep these values of Dm and R fixed, and vary #13. The Y 
matrix will change according to (3.1), but the Yukawa eigenvalues, yf will not change, as 
it is obvious from eq.(4.2). From the point of view of the Vl— parametrization the change 
in Y is due to a change in Vr and Vl- Therefore (Y^Y) 21 departs necessarily from its 
maximum value (4.6) and can only decrease in magnitude. However, from the point of view 
of the R— parametrization, the approximate expression eq.(4.3) tells us that an increase of 
S13 should reflect in an increase in the magnitude of (Y^Y) 21 ! So, at least we have found 
a choice of R for which the impact of #13 on (Y^Y) 21 [and thus of BR(/i — > e, 7)] goes 
exactly opposite that claimed in the literature. 

4.2 The reason behind 

The solution to the previous conflict can be found by doing the above steps explicit. We 
start with a choice for the Yukawa eigenvalues, Dy, and a Vl matrix that maximizes 
(Y^Y) 21 in eq.(4.5). The corresponding Dm, Vr can be obtained from eq.(3.4). One can 
construct now the Y matrix from (3.3). Then, using the R— parametrization (3.1), it is 
straightforward to derive the R— matrix (say R) that corresponds to this optimal choice: 

k = D ^M V ^ V L UD ^ > ( 47 ) 

(it is funny to check that R is orthogonal indeed) . If we keep now R and Dm constant but 
change the MNS matrix, U —> U' (e.g. by varying #13), it is straightforward from eq.(4.2) 
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that the new Y'^Y' matrix reads 

Y /1 Y' = U'U ] V L DlV ] L UU' ] = U'U ] Y ] YUU' ] . (4.8) 

Obviously, (Y'tY'^i < |(Y"^Y)2i| (recall that Vl was "designed" to maximize this quan- 
tity in eq.(4.5)). So something goes wrong with the argument used to obtain eqs. (4.3), 
(4.4). To see what, we construct the WDmR matrix that appears in the expansion (4.3): 

R^D M R = D^L UWYU D^L . (4.9) 

Hence 

1 



R ] D M R 



(X 



l 3 \[^i n d 

The consequence is that all terms neglected in (4.3) are in principle as large as the first 
term: the expansion is not sensible. In Fig. 1 we have plotted (|Y^Y)2i| 2 against #13, 
keeping R, Dm constant, for this particular example. As expected, the maximum occurs 
at #13 = 3° and for larger #13 the matrix element decreases. 



(4.10) 
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Figure 1. (|Y^Y)2i| 2 vs. #13 for a given R— matrix, correspondent the expression (4.7). 
Here D M = (10 10 , 10 11 , 10 12 ) GeV. 

Thus we have constructed an explicit counter-example, where the dependence of (Y^Y^i 
on #13 goes opposite that naively expected. This raises the question: Is this "wrong" be- 
haviour a consequence of the special choice of the seesaw parameters (R and Dm) taken 
above or it is more general? What can we expect for a generic choice of R and .Dm? 



5 Fair scans in the R— matrix and the perturbativity condi- 
tion 

The final questions of the previous section pose an interesting issue: how can we perform 
a truly generic scan in the R— matrix? In the literature it is common to separate the real 
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and imaginary parts of the Oi angles appearing in the parametrization (3.2). Then the real 
part is varied as an ordinary angle, say < Re Oi < 2tt, and the imaginary part is varied 
within a similar range. Note that such scan in Re Oi is really general, but the one in Im Oi 
is not. In principle Im Oi can take any value in the {—00,00} range. Of course a too-large 
value of Im Oi is not realistic, since it leads to non-perturbative Yukawa couplings, making 
the whole approach inconsistent. In this section we examine what restrictions does the 
perturbativity criterion impose on the magnitude of the Rij entries. We will see that they 
are quite simple, but very different from just choosing a certain range for Im Oi, which is 
the usual practice. 

The perturbativity requirement has to do with the Yukawa eigenvalues, Dy. Since, 
for a given Dm and R, these do not depend on the U— matrix [see eqs. (3.1) or (4.2)], the 
perturbativity criterion cannot depend on U either. A simple and sensible approach is to 
impose a constraint on the trace of Yukawa couplings, say 

trYtY = £y? <3 , (5.1) 

i 

(of course, any 0(1) number is as good as 3 here). Now, from eq.(4.2), 



trYtY = kArIDmR],. , (5.2) 

3=1,2,3 



.1,1 



so the perturbativity constraint (5.1) translates into 

\Rii\ 2 ~ tt— • (5-3) 

This condition is very handy and easy-to-use. Besides, it clearly applies whether or not 
we consider a supersymmetric version of the seesaw. An important remark is that the 
perturbativity requirement does not affect equally the magnitude of the various Rij entries. 
Actually, they can be easily different by orders of magnitude. This is clearly in contrast 
with typical scans of the R— matrix in the literature. It also explains the structure of the R- 
matrix [eq.(4.7)]. Recall that R was constructed to maximize (YtY) 21 at a certain value of 
#13. Then (YtY) 21 decreases for increasing #13, in contrast to the usual behaviour observed 
in the literature. But R is not a usual matrix considered in the literature. Actually, it 
"exploits" the perturbativity condition (5.3) to the extreme, as is clear from eq.(4.7). But 
it still corresponds to a perfectly sensible Y matrix. 

Now, we can pose the following question: For a given Dm, once R is scanned in all 
its generality (respecting perturbativity and orthogonality conditions), what is the corre- 
sponding range for (Y^Y) 21 and how does it change when #13 is varied? Intuitively, since 
yl < trYtY, we can expect a global range (see eq.(4.6)) 



0< 



YtY' 



21 



< 1 trY t Y . (5.4) 



But still we do not know how feasible or natural is to reach these bounds depending on 
the value of 6*13, or how (YtY) 21 changes with #13 for a fixed vanilla R. We have studied 
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this in a numerical way (more details in short), but we can get insight into these issues 
by examining the expression of (Y^Y) 21 in the R— parametrization [eq.(4.2)] more closely. 
This expression can be written as 

(YtY)^ = Y, U 2k ^R* jk Mj (iW^i* 3 + J2 R ^ U *i] ■ ( 5 - 5 ) 

j,k=l V i=l / 

Since \U\z\ = |si3 1 , the first term within the brackets is the responsible for the #13— dependence 
observed in the literature. However, the other two terms within the brackets can be easily 
as big as the first one. Note first that is only a factor ~ l/v6 smaller than a 
difference that is easily compensated by the fact that | 1 < \U\ 2 \. Besides, the pertur- 
bativity condition (5.3) implies that typically all the \J~k factors in (5.5) are compensated 
by the typical sizes of the R— matrix elements. In consequence, changing #13 is not likely 
to have a noticeable impact on (Y^Y) 21 , certainly not orders of magnitude for vanilla 
R— matrices. This simple argument can be made more rigorous (and cumbersome) once 
the orthogonality conditions on the R entries are imposed. But the basic result, that the 
range of (Y^Y) 21 cannot depend much on the value of #13, remains. This contradicts the 
common lore in the literature, and it is much more consistent with the result obtained 
using the Vl— parametrization. 

Let us now show the results of the numeric scan. First of all, we need a systematic pro- 
cedure to scan the whole range of R— matrices, consistent with the perturbativity condition 
(5.3) and the orthogonality condition, R T R = 1. A simple way to do it is explained in Ap- 
pendix A. Then, for each R— matrix considered, we scan #13 in the 0° — 10° range. Besides, 
for the numerical example we have taken the following values for the other parameters: 

m x = l(T 12 GeV, m 2 = 9 x 1(T 12 GeV, m 3 = 5 x KT 11 GeV, 
012 = 7r/6, 6*23 = vr/4, (5 = ^i = 2 = O, 

Mi = 10 10 GeV, M 2 = 10 11 GeV, M 3 = 10 12 GeV . (5.6) 

Furthermore, we have used tan/3 > 10, so that (H$) ~ v/V2, with v = 246 GeV. The 
results of the scan in R and #13 are shown in Fig. 2. As expected, the dependence of 
(Y^Y) 21 on #13 is very small, almost negligible. Note that, indeed, there are cases for 
which the dependence is stronger, corresponding to R— matrices which are far below the 
perturbativity limit (5.3), but those are statistically rare exceptions. 

6 The branching ratio BR(/j — > 6,7) 

In oder to translate these results about (Y^Y) 21 into predictions for the branching ratio 
BR(/j — > e, 7), we have to assume first a supersymmetric model. We have chosen a minimal 
supergravity (mSUGRA) model defined by 

m = 500 GeV, M 1/2 = 250 GeV, A = -100 GeV, tan/3 = 10, (6.1) 

where mo, M x / 2 and A$ are the universal scalar mass, gaugino mass and trilinear coupling 
at the unification scale Mx- Then we have calculated BR(/j — > e,7) for the same random 
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Figure 2. Scatter plot of (|Y^Y)2i| 2 vs. #13 for different R— matrices obeying the 
perturbativity condition (5.3). 

set of R— matrices used for Fig. 2. The computation was made by means of an own 
modified version of the SPheno code [8], which uses the full one-loop expressions of ref. [9]. 
The results are shown in Fig. 3 (left panel). As expected, the dependence of BR(/j — > e, 7) 
on #13 follows closely the one of |(Y^Y)2i| , shown in Fig. 2. Concerning the size of 
BR(/i —> e, 7), we see that in general is very large, quite above the experimental upper 
bound. This is in fact not surprising: using the approximate general range (5.4) and the 
approximate formula (4.1), one can check that the expected branching ratio is very large. 
Some comments are in order here. First, although it is common lore that BR(^ —> e, 7) can 
be quite large for a typical minimal SUGRA model, normally the actual values quoted in the 
literature are below those obtained here. This is because the scans in the R— matrix were 
not complete and many R— matrices compatible with perturbativity were not considered. 
Second, although it is not visible in Fig. 3, there are of course choices of R leading to 
branching ratios well below the experimental limit. What happens is that, scanning the 
space of the R— matrices in the way we did it, the number of those "good" R— matrices is 
relatively very small. We find this a very suggesting result. However, one has to keep in 
mind that scanning the R— parameter space in a different way (in Bayesian language, using 
a different prior for that space), the abundance of those "good" R— matrices will change. 
The dependence of these results on the prior is out of the scope of this paper. Nevertheless 
it is worth noticing that the same prior dependence was implicit in the scatter plots shown 
in the previous literature. Let us finally remark that by changing the parameters (6.1) 
of the mSUGRA model, the branching ratio changes parametrically, as indicated in the 
approximate formula (4.1). 

For the sake of comparison between parametrizations, we have also shown in Fig. 
3 (right panel) a similar survey using the Vl— parametrization. In this case we have to 
choose the values of the three neutrino Yukawa couplings. We have taken y\ = 0.0011, 
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yi = 0.03, 2/3 = 1. Besides, we have chosen the same mSUGRA model defined in eq.(6.1). 
Recall that in this parametrization (Y' Y) - does not depend at all on #13. The branching 
ratio does, but in a marginal way. All this is apparent from Fig. 3. Now, comparing the 
surveys with the R— and Vl— parametrizations (left and right panels) we note a similar 
insensitivity to #13, which is one of our main results. Besides, in the figure we see that 
the branching ratios in the R— parametrization can be more than one order of magnitude 
larger than in the Vl one. This is mainly due to the choice y% = 1 in the latter. Note that 
in the R— parametrization we have imposed trY^Y < 3, which allows ?/| < 3. Therefore 
the upper limit of |(Y^Y)2i| 2 can be almost one order of magnitude bigger than in this Vl 
survey, see equations (5.4), (5.4). Furthermore, in the Vl— parametrization the righthanded 
neutrino masses are an output. In the example chosen they come out typically bigger (and 
less degenerate) than the choice made for the R— parametrization, eq.(5.6). This pushes 
downwards further the branching ratio through the log factor, see eq.(4.1). Hence, the two 
surveys are perfectly consistent. 




02468 10 02468 10 



o 13 [grad] »1 3 farad] 

Figure 3. BR(fi — > e, 7) vs. #13. (left) R— parametrization; (right) Vl— parametrization. 



Of course, the results of this section could be different if the set of R— matrices (or Vl 
and Vr matrices for the V^-parametrization) is constrained by additional considerations. 
This is the case of GUT models (as those studied in [10]), where non-trivial correlations 
can indeed occur. 

The last point raises a final question: what happens if leptogenesis constraints are 
imposed in the scenario? Actually, the impact of #13 on LFV processes was reported to get 
reinforced once successful leptogenesis is incorporated in the analysis. In the next section, 
we re- visit the leptogenesis issue. 



7 Constraints from leptogenesis 

The baryon asymmetry of the universe (BAU) is usually defined as the ratio of the number 
density of baryons to the number density of photons re 7 [11]. Its present experimental 
value is [12] 

— = (6.19 ±0.15) x KT 10 . (7.1) 
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Perhaps the most popular mechanism to generate such BAU is nowadays thermal lepto- 
genesis, in which a net lepton number is produced by the out-of-equilibrium decay of the 
(seesaw) right-handed neutrinos. Then the lepton number is converted into baryon num- 
ber by sphaleron-mediated processes. Note here that the relevant Lagrangian, defined by 
the superpotential (2.1), contains the required lepton number and CP violating-terms (the 
latter are provided by appropriate phases in the Y matrix). 

The final value for the BAU is given by: 



n B 



1U 



Csphal e 1] 



(7.2) 



where n\ is the equilibrium number density of the lightest righthanded neutrino, which 
is the main responsible for the asymmetry for hierarchical righthanded masses, M\ -C 
M2 -C M3; C sp hai contains the sphaleron effect, e is the CP- violating contribution to the 
asymmetry, and 77 is the efficiency factor, which takes into account the partial erasure of the 
CP asymmetry by inverse decays and scattering processes. In the minimal supersymmetric 
standard model eq.(7.2) can be written as [11] 



n B 



-1.04 x 1CT 2 e 77 . 



(7.3) 



The thermal production of righthanded neutrinos is suppressed unless M\ < Tr, where 
Tr is the reheating temperature after inflation. On the other hand, Tr cannot be much 
larger than 10 10 GeV, to avoid the gravitino problem. In the following we will assume 
Mi ~ Tr ~ 10 10 GeV. In this temperature regime, one or more charged-lepton mass- 
eigenstates £ (£ = e, /j, r) are in equilibrium in the thermal bath, and flavour effects can be 
significant because the corresponding lepton asymmetries follow an independent evolution 
[13-15]. However we will not consider flavour effects here, since our main goal in this paper 
is to examine the dependence of the results on #13, and to compare the results with the 
previous literature (where flavor effects were not considered) 1 . 

In the unflavored case, the efficiency factor, 77, is given by 



2(m e + m M + m T ) 



???* 



1.16 



+ 



2m* 



rh e + rh u + rn T 



(7.4) 



where 



rhi = |Yii| 2 ^/ Ml 

m* = 4vr7j 2 Fi/M 1 2 
= H(T = Ma) = 1.66V^M 2 /M Planck , 



1 = e, /j, r 
l(T 3 ey 



A^Imssm = 228.75 



(7.5) 
(7.6) 



1 It is worth mentioning here that in [16] it was shown that successful lcptogenesis is possible 
within the SUSY seesaw for any value of the still unmeasured low energy neutrino parameters 
(including #13), taking into account flavour effects and using the Vl - parametrization. This may 
be an indication that flavour effects are not going to introduce any dramatic dependence on #13. 
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with H denoting the Hubble parameter. The CP factor, e, is given by [17]: 



1 E^rMK^JuH 9(Mf/M?) 



8tt 



(YYt) 



11 



(7.7) 



where g(x) 



In 



x + 1 



x 



Note that both r\ and e depend on the Yukawa couplings through the combination 
YYt, which, in the R— parametrization -eq.(3.1)- is independent of U 



YYt = D ^ m RD K R^D 



M 



(7.8) 



Hence, the BAU, given by (7.3), does not depend on #13. 

Still, it could happen that leptogenesis constraints strengthened the dependence of 
BR(/x —7- e,7) on #13. As we have seen in the previous section, there are some R— matrices 
for which the dependence of (YtY) 21 (and thus of the branching ratio) on #13 is impor- 
tant. If (in an extreme case) those R— matrices were precisely those selected by successful 
leptogenesis, we would find a strong dependence of the branching ratio on #13 in the com- 
plete scenario. To analyze whether this possibility (or a more moderate one) takes really 
place, it is important to perform a complete scan of the R— matrix space. Performing 
partial scans in this space can be useful to show particular features, but it may introduce 
unwanted biasses: one could artificially select R— matrices that lead to strong dependences 
of BR(/i — > e, 7) on #13 (as has happened in previous literature). 




1e-12 1 1 1 1 1 1 

2 4 6 8 10 

013 [g rad ] 



Figure 4. BR(/j — > e, 7) vs. #13 in R— parametrization; including leptogenesis constraint. 

Fig. 4 shows the dependence of the BR(/i — > e, 7) on #13 after imposing successful lep- 
togenesis and scanning R in its whole parameter space. Clearly, no important dependence 
on #13 is observed. The trend is very similar to the one with no leptogenesis constraints 
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(Fig. 3). The conclusion is that leptogenesis constraints do not enhance (or create) any 
dependence on #13. 

Note that when the leptogenesis constraint is imposed, the values BR(/j —> e, 7) de- 
crease at least one order of magnitude with respect to the case without leptogenesis, Fig. 3 
(left). This behaviour can be understood from eqs. (7.3), (7.7) and (7.4). The baryon 
asymmetry depends on the neutrino Yukawa matrix Y both through the CP asymmetry e 
and the efficiency factor rj. Note from eq. (7.7) that the CP asymmetry is mainly driven 
by the ur 2 and vr s Yukawa couplings (that is, the second and third rows of R), since the 
dependence on the vri Yukawas approximately cancels out due to the (YY')n factor in 
the denominator. As a consequence, we find a mild dependence of e on the elements in the 
first row of R (see Fig. 4). 

The efficiency factor, rj, in eq. (7.4) smoothly interpolates between the weak rhi <C 
m*) and strong i^2 i rhi 3> m*) washout regimes. However for all the points in our scans 
we find that Y2i ™i m * > thus the efficiency factor is given by 



V 



2T m l" 1 - 16 

' , <• <.I<-T. (7.9) 



and it becomes clear that 77 decreases for larger Yij, and thus for larger R— matrix elements, 
see eq. (3.1). So, after imposing enough BAU, only those R— matrices which are sufficiently 
small in order to optimize rj are selected, thus reducing the final value of the branching 
ratio. 




Figure 5. Scatter plot of e and 77 for random R-matrices, every point represents a 
different R— matrix. The color coordinate represents the average absolute-value of the 
elements in the first row of R. The red and green lines denote the lower and upper BAU 
limits, respectively. For clarity in the plot, we have relaxed the allowed BAU window to 
be [5 - 7] x 1(T 10 . 
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In Fig. 5 we have performed a scatter plot showing the values of r/ and e. The exper- 
imental window is also shown for reference. From the figure it is clear that enough BAU 
is only produced if the elements of the first row in R are smaller than their perturbative 
bound, in agreement with our previous analytical estimations. This is the reason for the 
decrease of the branching ratio. Besides, there is no color gradient in the e— direction, 
which means that the same values of e can be reached for any value of the first row entries 
of the R— matrix. This is in contrast with the case of rj, which changes dramatically when 
varying the typical size of the first row R— matrix elements. 

8 Conclusions 

The main results and conclusions of this paper are the following: 

• The extended believe that the branching ratio BR(/x — > e, 7) in supersymmetric 
seesaw models depends strongly on the value of #13 does not hold after a careful 
analytical and numerical study. 

• We have analyzed this issue using two alternative parametrizations of the 9 degrees 
of freedom that, besides the 9 low-energy observables (neutrino masses, mixings and 
phases), expand the parameter space of the seesaw scenario. This amounts to two 
alternative ways of traveling across this 9-dimensional space or, in other words, of 
parametrizing our ignorance. These are called the i?-parametrization (sect. 3.1) and 
the Vk-parametrization (sect. 3.2) 

• The main potential dependence of BR(/j — > e, 7) on #13 occurs through the Y^Y^ ma- 
trix, where Yj, is the neutrino Yukawa matrix. In the V^-parametrization, this quan- 
tity is trivially insensitive to #13 (or to any observable parameter), so BR(/i — > e,j) 
is. In the -R-parametrization (yJ/Y^)^- has a dependence on #13, which essentially 
disappears once the 9-dimensional parameter space is fairly covered. 

• In the ii-parametrization (which is the most common in the literature) a fair scan 
implies to allow all the complex i?-matrices compatible with orthogonality and per- 
turbativity of the Yukawa couplings. The latter requirement has not been properly 
taken into account in former literature. As a consequence previous scans in the 
space of the i?-matrices are typically biassed, since many possible i?-matrices were 
excluded from the beginning. 

• We give a very simple rule to incorporate the perturbativity of Yukawa couplings as a 
condition in the entries of the orthogonal i?-matrix. It is given by eq. (5.3). We also 
give (in Appendix A) an straightforward procedure to completetly scan the space of 
complex i?-matrices in a consistent way with this requirement and the orthogonality 
one. 

• Once such scan is performed the branching ratio BR(/i — > e, 7) gets very insensitive to 
#13, as already mentioned. Moreover, the values of the branching ratio are typically 
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larger than those quoted in the literature. This comes from the fact that many 
possibilities that were disregarded from the beginning (typically with large i?-entries 
and thus sizable Yukawas) turn out to be perfectly compatible with the perturbativity 
constraint, which had not been taken into account. 

• We find this increase of BR(/i — > e, 7) a very suggesting result. However, one has to 
keep in mind that if we scanned the R— parameter space in a different way (though 
still covering the whole space) the relative abundance of points with large BR(/x — > 
e, 7) might change. In Bayesian language, this is equivalent to use a different prior 
for the parameter space. The dependence of the typical size of BR(/i — > e, 7) on the 
prior is out of the scope of this paper. Nevertheless it is worth noticing that the 
same prior dependence was implicit in the scatter plots shown in previous literature. 

• Including leptogenesis constraints (disregarding flavour effects) in the analysis does 
not introduce any further dependence of BR(/i — > e,j) on #13. The main impact of 
leptogenesis, besides remarkably reducing the acceptable volume of the i?-parameter- 
space, is a decrease of BR(/x — > e, 7) by more than one order of magnitude. This 
comes from the fact that the efficiency factor, rj, decreases for large .R-matrix entries. 
Hence, successful leptogenesis prefers smaller ones, and thus smaller BR(/i — > e, 7). 

As a concluding remark, scanning the parameter space of the R— matrix in the full allowed 
range is necessary in order to make general statements about predictions of the seesaw 
scenario. This has not been taken into account in former works, and we consider it an 
interesting line of work to be explored in forthcoming projects. 
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A General scan of the R— matrix 



In this section we explain the details of the scan made on the R— matrices, to cover all 
possibilities compatible with orthogonality, R T R = 1, and the perturbativity condition 
(5.1, 5.3), which for convenience we repeat here: 

The algorithm presented below has been designed to be easily modified if one considers 
also leptogenesis constraints. 

In general, given a normal hierarchy among neutrinos, K\ <C k<i <C K3, condition (A.l) 
tells us that the R— matrix elements are allowed to be larger (in absolute value) when you 
move from bottom-right to top-left in the matrix. Thus, the element .R33 (-R11) presents the 
smallest (largest) upper bound. On the other hand, orthogonality implies that, in practice, 
not all the elements of R— matrix can reach their corresponding perturbativity limit. Nor- 
mally, if one sets any entry of the first row (column), at its maximum magnitude, then the 
corresponding entries in the same column (row) cannot satisfy orthogonality without violat- 
ing the perturbativity bound (A.l). Conversely, if the matrix elements R22, R23, R32, R33, 
satisfy their perturbativity constraints, then normally the entries of the first row and col- 
umn would do it as well. Hence, it makes sense to start imposing perturbativity in the 
the bottom-right part of the R— matrix. Then the rest is constructed automatically from 
orthogonality, requiring a final cross-check of the perturbativity condition. 

More in detail, using the parametrization (3.2), we see that 

-^32 + -^33 = C 2 > 

but the perturbativity condition tells us that R\ 2 is allowed to be much larger than -R332, 
up to phases. So, we can use the perturbativity upper bound on .R32 to scan C2 

c 2 = |c 2 |e^ 2 , |c 2 | < j= , 4>2 G [0,2tt]. (A.2) 

V M3K2 

Here the phase <f>2 (and those appearing below) is assumed to be a random number within 
its interval. Now, for each value of C2 we derive S2 (obviously, up to the sign) and exploit 
the perturbativity condition on i?33 to scan c\\ 

Cl = | Cl |e^, |d| < <£i<E[0,27r] (A.3) 

\c 2 \VM 3 k 3 

(again, for each value of c\ there are two values of si). At this point we have used the 
bounds on -^32,-^33, and the orthogonality condition, to scan #1,6*2. 

Now, we have to scan the third complex angle, #3. A convenient way to do it is by 
using the -R23 element (recall we prefer to impose perturbativity in the bottom-right part 
of R). Thus we scan 

R23 = |«2 3 |e^ 23 , |«2 3 | < -7=t=, 023 G [0,2vr] . (A.4) 

V M2K3 
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For each value of R23 we derive the two possible values of S3: 



_ S 2 ^23Cl ± y/sfj-R^ + 4 + cjsj) 

33 ~ 44 + 4 ' { ' 

Once again, for each value of S3 there are two of C3. Finally, we cross-check for each of 
these values that the corresponding R— matrix is indeed consistent with the perturbativity 
condition. One can do that directly by using equation (5.1) or by checking eq. (A.l) for 
the remaining entries. 

When leptogenesis constraints are included, most of the initially-allowed values for the 
first-row entries of R, become too big, since they lead to a small efficiency factor, r] (see 
sec. 7). In this case, it pays to constrain from the beginning any of the entries in the first 
row, e.g. R13. Then, instead of the previous scan in R23, one scans -R13 as: 

Rl3 = l^isle^ 13 , |#13|<1, 013 6 [0,2tt]. (A.6) 

Note that this upper bound on | R13 | is normally much lower than its perturbativity bound 
(~ 10 2 in our case), but still is way larger than required by leptogenesis constraints. Now, 
for each value of R13 the value of S3 is given by 

R 13Sl ±y/c(s 2 2 (-R( 3 + 4+cisf) 

33 = 4Tc&2 ' (AJ) 

which replaces eq.(A.5). 

As a final comment, recall that the parametrization of R given by eq. (3.2) is completely 
general up to reflections changing the sign of det R. In our case, however, the scan is 
completely general since all the relevant physical quantities are invariant under global 
changes of the sign of R. 



B Details of the numerical computations 

B.l R— parametrization 

In this section we explain the strategy for computing the branching ratio BR(/j — > e, 7) for 
each point in the scan of the R— matrices (see Appendix A). 

We have adopted an mSUGRA framework, with universal soft terms at the GUT scale, 

M X , 

(m£)ij = (m^Jjj = m l, (A e )»j = A (Y e ) ij , (B.l) 

where m^, rrig R and A e are the left- and right-handed slepton mass-squared matrices, and 
the matrix of slepton trilinear couplings. At Mx all the soft terms are diagonal in the 
basis in which Y e is diagonal. Below Mx, the RG running of produces off-diagonal 
entries, due, essentially, to the contribution proportional to Y^Y V . This constitutes the 
main source of flavor violation. 
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We work with our own modified version of the SPheno [8] code. The matrix is of 
special interest since, among the slepton mass matrices it is by far the one that develops 
larger off-diagonal terms. At the seesaw scale, M, we evaluate as 

ml = D m£ - ^(3mg + 4j)YtD L Y„ . (B.2) 

Here D m 2 is the result of running m£ from Mx until M, switching off the contribution 
from neutrino Yukawa couplings. The second term in (B.2) is the contribution coming 
from the neutrino Yukawas, evaluated at the leading-log approximation. This contribution 
contains the off-diagonal entries of m^. The value of Y y at the M— scale is obtained from 
the R— parametrization formula (3.1) for each point in the scan of the R— matrices. Finally, 
is run down to low-energy (neutrino Yukawas do not play any role in this RG-interval 
since right-handed neutrinos are decoupled). 

The rest of physical quantities (charged slepton mass matrices, gauge couplings, GUT 
scale, charged lepton yukawas, etc.) are taken directly from SPheno, which imposes the 
Mx scale to be the one where gauge couplings unify. We also extracted from SPheno the 
parameters of the neutralinos and charginos. Finally, we followed ref. [3] (implemented in 
SPheno) to calculate the branching ratio at 1-loop level. 

B.2 Vl— parametrization 

For the Vl— parametrization, we have used the SPheno code as well. However, the original 
code is not prepared to introduce the initial parameters according to this parametrization. 
In particular, the original version works with given values of right-handed neutrinos masses, 
Mj, from the beginning. But in the Vl— parametrization, Mj (and Vr) are obtained at the 
seesaw scale, M (a suitable average of Mj), from eq. (3.4). On the other hand, the neutrino 
Yukawa eigenvalues, Dy and the Vl matrix are indeed initial parameters, defined at the 
high-scale (and these are the ones in which we perform our scan of the parameter space). 

Consequently, we modified the code, incorporating an iterative procedure to determine 
M{ in a consistent way with all the boundary conditions (at low- and high-scale). In 
this way, the full matrix is obtained directly from SPheno, and can be used to the 
computation of BR(/i — > e, 7). 
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